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Several  missiles  and  projectiles  have  periods  during  their  flight  where  the  body  is  spinning  and 
must  be  controlled.  It  has  been  stated  that  a  roll  stabilized  autopilot  may  be  more  appropriate 
during  this  phase  of  the  flight.  This  paper  first  examines  the  full  nonlinear  equations  of  motion  in 
both  the  roll  stabilized  frame  and  the  body  axis  system  assuming  a  rigid  axially  symmetric  body. 
Two  and  three  loop  autopilots  are  developed  using  an  optimal  control  approach  for  a  non  rolling 
airframe.  These  autopilots  are  examined  to  determine  how  well  they  are  able  to  follow  inertial 
commands.  An  extremely  simplified  example  is  presented  so  that  the  basic  characteristics  of  the 
closed  loop  system  can  be  examined.  The  results  are  presented  and  conclusions  drawn.  Effects 
of  different  acceleration  levels,  including  the  axial  moment  of  inertia  and  very  high  roll  rates  are 
examined. 


I.  Introduction 

Several  missiles  and  projectiles  have  periods  during  their  flight  where  the  body  is  spinning  and  must  be 
controlled.  This  paper  develops  the  full  nonlinear  equations  of  motion  in  both  the  body  and  what  is  called 
the  roll  stabilized  coordinate  frames.  The  roll  stabilized  frame  is  a  frame  where  the  Euler  roll  rate  is  always 
zero.  Four  autopilot  topologies  were  examined  for  control  of  a  rolling  airframe.  These  were  the  classical  three 
loop  implemented  in  the  body  axis  system,  the  classical  three  loop  implemented  in  a  roll  stabilized  system 
and  a  two  loop  autopilot  implemented  in  both  the  body  and  roll  stabilized  frames.  It  is  shown  the  two  loop 
implementations  are  identical.  There  is  no  advantage  of  using  the  roll  stabilized  coordinate  frame  for  the 
two  loop  design.  The  three  loop  autopilots  produce  significantly  different  results.  An  extremely  simplified 
numerical  example  is  presented  so  that  the  basic  characteristics  of  the  closed  loop  system  can  be  examined. 
Results  are  presented  and  discussed.  Finally  conclusions  are  reached. 

II.  Coordinate  Frames  and  Notation 

There  are  three  coordinate  frames  used  in  the  derivation  of  the  equations  of  motion.  There  is  the  inertial 
frame  designated  with  and  the  roll  stabilized  frame  labeled  “</>”  and  the  body  frame  designated  UB” . 
The  inertial  frame  is  the  frame  in  which  Newtons  Equations  of  Motion  can  be  applied.  The  roll  stabilized 
frame  is  defined  such  that  the  Euler  angular  rate  is  zero  for  all  time  (4>  =  0).  The  body  frame  is  attached 
to  the  spinning  body.  The  unit  vector  triplets  (ix,jx,  kx)  will  be  used  for  each  coordinate  system.  The 
rotation  sequence  from  the  inertial  to  the  roll  stabilized  frame  is  first  through  a  rotation  about  the  inertial 
z  axis  through  the  heading  angle  T.  Using  an  intermediate  coordinate  system  “1”.  The  second  rotation  is 
about  the  intermediate  y  axis  through  the  angle  0.  The  direction  cosine  matrix  is 
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The  angular  velocity  between  the  systems  (read  as  the  angular  velocity  of  frame  “0”  relative  to  “7” )  is 

Z+n  =  Tfc/  +  0j0  =  pi^  +  qjj,  +  rk 4 

The  body  axis  system  is  arrived  through  one  more  rotation  about  the  roll  stabilized  x  axis.  That  is 


lB 

^<j) 

JB 

— 

3<t> 

k b 

k<f, 

where 


C 


B/4> 


0 

cos(3>) 


0 

sin(<f>) 


0  —  sin($)  cos(3>) 

and  the  angular  velocity  of  the  body  frame  relative  to  the  roll  stabilized  frame  is 

^B/c/l  =  <f>*<p 

There  is  absolutely  nothing  special  about  the  above  definitions.  The  exact  same  rotation  sequence  and 
equations  can  be  found  in  Roskarn,  Blakelock  and  Etkin.  The  only  difference  is  the  roll  stabilized  frame  was 
defined  as  an  intermediate  frame.  The  angular  velocity  of  the  body  frame  relative  to  the  inertial  frame  is 

UJb/I  =  <^B/c/>  +  &4>/I  =  ^kj  +  +  &lB  =  PlB  +  Q]b  +  RkB 

III.  Linear  Momentum  and  the  Force  Equation 

The  linear  momentum  of  a  rolling  rigid  airframe  is 

P  =  mIVC/o 

where  the  velocity  is  defined  as  the  time  rate  of  change  of  the  position  of  the  center  of  mass  “C”  relative  to 
an  inertially  fixed  point  “O”  as  seen  by  an  inertial  observer.  The  notion  will  be 


vc/o  = 


IdRl 


c/o 


d  t 


Newtons  Law  states  that  the  mass  times  the  acceleration  is  equal  to  the  sum  of  the  applied  external  forces. 
That  is 

'd’Vc/o 


FExt  = 


ni¬ 


di 


The  derivative  can  be  expanded  using  the  transport  theorem,  or 

? c/o  _  BdIVc/o 


7d  IVn 


+  &B/I  x  Vc/O 


dt  dt 

IV.  Angular  Momentum  and  the  Moment  Equation 

The  angular  momentum  about  the  airframe  center  of  mass  of  a  rotating  rigid  airframe  is 

H  =  IC  ■  Cjb/i 

For  an  axially  symmetric  body,  the  moment  of  inertia  tensor  reduces  to 

Ic  =  IxxlBlB  +  IjBjB  +  IkBks 

This  notation  may  look  a  little  strange,  it  is  simply  a  tensor  notation,  the  results  are  the  same  as  for  the 
standard  matrix  approach  only  a  matrix  cross  product  does  not  need  to  be  defined.  For  a  rigid  body  Newtons 
law  reduces  to  the  change  in  angular  momentum  is  equal  to  the  sum  of  the  external  moments,  or 

'dH 


d( 
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again  using  the  transport  theorem 


'AH  BdH  w 

+  UJ,3/r  x  H 


V.  Body  Axis  System  Equations  of  Motion 

The  equations  of  motion  in  the  body  axis  system  are  derived  using  the  body  axis  variables,  the  linear 
velocity  is 

'Vc/o  =  UiB  +  VjB  +  WkB 

and  the  angular  velocity  is 


&B/I  =  PiB  +  Qjb  +  Rks 


The  force  equation  is 


=  Uib  +  Vjs  +  WkB  +  (PiB  +  Qjb  +  Pkis^j  x  (uiB  +  VjB  +  Wk/^j 


or  in  component  form 

U  +  QW-RV  =  — 

m 

V  +  RU  —  PW  =  — 

m 

W  +  PV-QU=  — 
m 

where  the  force  components  are  in  the  body  axis  system. 

The  moment  equation  is 


=  ~aT 


—  {j-xxP  +  IQR  —  IQ  IIs)  %b  +  {jQ  +  IxxP R  —  IP r)  jb 
+  ^ IR  +  IPQ  —  IxxPQ^j  kB 


or  in  component  form 


IxxP  =  Mx 

IQ  +  IxxP  R  —  IP  R  =  My 


IR  +  IPQ  —  IxxPQ  —  Mz 

where  the  moment  components  are  in  the  body  axis  system. 


VI.  Roll  Stabilized  Equations  of  Motion 

The  equations  of  motion  in  the  roll  stabilized  axis  system  are  derived  using  the  roll  stabilized  axis 
variables.  The  linear  velocity  is 

rVc/o  =  +  vj $  +  wk $ 

and  the  angular  velocity  is 


+  qjc/,  +  rkfj, 


The  force  equation  is 


^  ^xt  =  +  vjtj,  +  wktf,  +  (pi^  +  qjtj,  +  rk ^  x  (ui^  +  vj $  +  wk </, 
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or  in  component  form 


u  +  qw  —  rv  = 


fx 

m 


v  +  ru  —  pw  = 


fy 

m 


w  +pv  —  qu  = 


fz 

m 


where  the  force  components  are  in  the  roll  stabilized  system.  Clearly  the  force  equation  have  not  been 
simplified  using  the  local  level  frame,  though  it  should  be  noted  the  roll  rate  in  the  roll  stabilized  system  is 
zero  at  zero  pitch  attitude  and  is  significantly  smaller  for  rolling  airframes  compared  to  the  roll  rate  of  the 
body  axis  system.  One  inconvenience  is  that  the  forces  are  usually  written  in  the  body  axis  system  but  for 
a  symmetric  body  this  difficulty  should  be  able  to  be  overcome. 

The  moment  equations  are 


MExt 


'dH 
d  t 


i>dH  _ 

IT  +  ^  x 


H 


*dl, 


C  •  Ub/I 

d  t 


+  tip/!  X  (ic  • 


we  must  now  pause  and  look  at  the  time  rate  of  change  of  the  inertia  tensor  as  seen  by  someone  in  the  roll 
stabilized  frame.  Fortunately  the  body  is  axially  symmetric  thus  using  the  coordinate  transformation  of  a 
tensor  we  find  the  moment  of  inertia  tensor  is 


Ic  —  I xx.i  <p  i ft  I  I 3  P3P  I  Ikpkp 

so  that  it  is  also  invariant  seen  by  the  roll  stabilized  frame,  thus 


Ext 


=  Ic 

=  Ic 


i’d  £3 


B/I 


d  t 

'Bdu 


d  t 


+  tip/I  x  (lc  ■  Ub/i') 

- - +■  tiB/p  x  Ub/i  \  +  tip/ 1  x  (Ic  •  tiB/i^ 


now  we  must  either  put  the  local  level  angular  velocity  or  the  body  angular  velocity  vector  in  a  different 
frame  so  that  the  cross  products  can  be  taken.  Choosing  the  later 

tiB/i  =  PiB  +  QjB  +  RkB 

=  Pip  +  [Q  cos(d>)  -  R  sin  ($)]  jp  +  [Q  sin($)  +  R  cos  ($)]  kp 

therefore 

^  '  MExt  —  (ixxipip  T  I/jp3p  T  Ikpkp'j  •  (PzB  T  Q3b  I  Rks'j 
T  (ixxipip  T  Ijp3p  T  Ikpkp)  *  *&ip  x  ujb/i 
+  (pip  +  Q3p  +  rkp'j  x  (lxxipip  +  Ijpjp  +  Ikpkp'j  ■  tiB/i 
Substituting  using  the  following  identities 


q  =  [Q  cos(d>)  —  R  sin  ($)] 
r  =  [Q  sin(<l>)  +  R  cos  ($)] 
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q  =  Qcos(3>)  —  Psin  ($)  —  TQ  sin(<I>)  —  <f>P  cos  ($) 


and 


r  =  Q  sin(<I>)  +  R  cos  (3>)  +  <1 )Q  cos(3>)  —  4>P  sin  ($) 
yields  the  component  equations 

I :i:x  V  =  'W'x 


Iq  +  IxxPr  -  Ipr  =  my 
If  +  Ipq  -  IxxPq  =  mz 


where  the  moments  are  written  in  the  roll  stabilized  frame. 

Notice  most  of  the  cross  coupling  is  gone  but  not  all  the  roll  rate  of  the  roll  stabilized  frame  (small  but 
usually  nonzero)  still  produces  a  cross  coupled  term  and  the  gyroscopic  effects  still  couple  the  pitch  and  yaw 
channels.  The  fact  that  the  roll  rate  of  the  roll  stabilized  frame  (p)  is  significantly  smaller  then  the  roll  rate 
of  the  body  axis  system  (P)  is  the  reason  to  hope  for  better  closed  loop  performance  using  the  roll  stabilized 
frame.  The  actual  simulations  that  will  be  used  to  produce  the  results  presented  later  always  use  the  body 
axis  system  of  equations  because  the  fins  and  aerodynamics  are  truly  attached  to  the  body. 


VII.  Alternate  Autopilots 

Four  autopilots  will  be  discussed.  They  are:  a  three  loop  in  the  body  axis  system;  a  three  loop  in  the  roll 
stabilized  frame;  a  two  loop  in  the  body  axis;  and  a  two  loop  in  the  roll  stabilized  frame.  Each  will  discussed 
in  detail. 


A.  Three  Loop  in  the  Body  Axis  System 

The  equations  for  the  fin  commands  for  the  three  loop  autopilot  using  body  axis  variables  are: 

5y  =  KIAy  J  ( AVm  -  KssAVc)  di  +  K,lf!  J  Rm  dt  +  KrRm 

fip  =  P-IAz  J  (KssAZc—AZm)  dt  +  Kg  J  Qm^t  +  KqQm 

The  subscript  “rn”  is  for  measured.  That  is,  these  are  measured  quantities.  This  paper  assumes  a  perfect 
Inertial  Measurement  Unit  (IMU)  in  that  the  IMU  is  at  the  center  of  mass  so  there  are  no  moment  arm  and 
there  are  no  dynamics  associated  with  the  measurement  device.  It  should  be  noted  the  output  of  the  IMU 
is  G’s  for  acceleration.  The  block  diagram  of  the  pitch  channel  is  presented  in  Figure  1. 


Figure  1.  Three  Loop  Topology 
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B.  Three  Loop  in  the  Roll  Stabilized  System 

The  equations  for  determining  commands  for  this  autopilot  are  first  the  deflections  in  the  roll  stabilized 
system  are  computed  from 

Sy  RlAy  j  (^Aym  ^ss^-y^  dt  d-  ^  V m  df  d~  RrVm 

Sp  =  RlAz  J  (^R-ssAZc  dt  T  I\.Q  J*  Qm  dt  RqC[m 

Then  these  are  transformed  into  the  body  fin  deflections  through 

5 Y  =  SY  cos (</>)  +  Sp  sin(^) 

5p  =  —Sy  sin((/>)  +  Sp  cos((f>) 

It  should  be  noted  the  particular  sign  convention  is  used  because  the  fin  deflections  are  not  a  vector  but  in 
order  to  get  the  correct  force  or  moment  out  of  the  fin  the  above  signs  must  be  used.  The  roll  stabilized 
measurement  variables  are 

AVm  =  AVm  cos  {(j))  +  AZm  sin(^) 

AZm  =  -AVm  sin(</>)  +  AZm  cos ((f>) 

and 

qm  =  Qm  cos{(/))  +  Rm  sin (</>) 
rm  =  -Qm  sin(^)  +  Rm  cos(^) 

The  whole  idea  of  this  autopilot  is  to  compute  the  fin  deflections  in  the  roll  stabilized  frame  and  modulate 
the  output  into  the  missile.  The  block  diagram  is  the  same  as  for  the  body  axis. 

C.  Two  Loop  in  the  Body  Axis  System 

The  equations  for  the  fin  commands  for  the  two  loop  autopilot  are: 

Sy  =  R Ay  ( Aym  KssAyc)  d“  Rj-qRrn 

Sp  =  Kaz  ( KssAZc  —  AZrn  )  +  KqQm 
The  block  diagram  of  the  pitch  channel  is  presented  in  Figure  2. 


Figure  2.  Two  Loop  Autopilot  Topology 
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D.  Two  Loop  in  the  Roll  Stabilized  System 

The  equations  for  the  fin  commands  for  the  roll  stabilized  two  loop  autopilot  are: 

^ Y  =  R  Ay  ^Aym  A ssAy\  -f-  Rrvm 

5p  =  Kaz  (kssAZc  —  AZrr\  +  Kqqm 

Sy  =  Sy  cos (</>)  +  5p  sin(0) 

Sp  =  —Sy  sin(</>)  +  Sp  cos((/>) 

again 

AVm  =  AVm  cos(</>)  +  AZm  sin ((/>) 

AZm  =  —AVm  sin  (cf>)  +  AZm  cos  (</>) 

and 

qm  =  Qm  cos{(/))  +  Rm  sin(^) 
rm  =  ~Qm  sin {</))  +  Rm  cos (0) 

The  block  diagram  is  the  same  as  for  the  body  axis. 


VIII.  Simplified  Example 


The  following  simplifications  were  made  in  order  to  understand  the  basic  dynamics  of  the  closed  loop 
systems.  The  effects  of  gravity  are  ignored,  the  x  axis  moment  of  inertia  is  zero  ( Ixx  =  0).  The  forward 
velocity  is  constant,  the  pitch  and  heading  angles  are  ignored  and  initially  zero,  the  roll  rate  is  constant. 
The  y  and  z  force  was  a  linear  function  in  sideslip  or  angle  of  attack  respectively,  The  y  and  z  moments  are 
linear  functions  of  both  angle  of  attack  or  sideslip  and  pitch  and  yaw  fin  deflections.  The  dynamics  were 
always  computed  in  the  body  axis  system.  The  dynamic  equations  are: 

a  =  tan  1  —  1  pM  =  tan 

Fx  =  0  Fy  =  YpMfiM  Fz  =  Zaa 
Mx  =  0  My  =  Maa  +  Mgp5p  Mz  =  NpMj3M  +  Nsy5y 


U  =  0 


P  =  0 


V=FA 

m 

qPA 

Ax  =  — 

mg 


RU  +  PW 


W  =  — 
m 


PV  +  QU 


M. 

PR  R  =  — — ^  —  PQ 


=  P 


Ay  =  ^L 
mg 


mg 


Equations  1  and  2  are  integrated  in  the  nonlinear  dynamics  model  to  produce  the  states. 
The  aerodynamics  are  further  simplified  by 


(1) 

(2) 


YpM  —  Fa  Ma  —  N/3m  MsP  —  NgY 
The  specifics  of  the  dynamics  model  are 


Variable 

Value 

Description 

Value 

m 

46.0 

Mass 

Kg 

i 

22 

Moment  of  Inertia 

Kg-m2 

Za 

-4712.4 

Force  due  to  Angle  of  Attack 

N 

Ma 

-549.8 

Moment  due  to  Angle  of  Attack 

N-m 

MsP 

-1570.8 

Moment  due  to  Fin  Deflection 

N-m 

U 

500 

Forward  Velocity 

m/s 

g 

9.81 

Gravity  Constant  used  for  Unit  Conversion 

m/s2 
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Roll  Stabilized  Acceleration 


Figure  3.  Roll  Stabilized  Acceleration 


The  autopilots  were  designed  using  a  linear  optimal  control  method.  The  specific  gains  for  the  three  loop 
designs  are 


KIAy  =  1.3671  Kg  =  2.4020  Kq  =  0.2594  Kss  =  1.0345 
and  for  the  two  loop  designs  are 

K  i  Ay  =  0.1650  I\q  =  0.2197  Kss  =  1-2293 

IX.  Results 

The  first  check  was  to  make  sure  the  autopilot  design  was  implemented  properly  in  the  nonlinear  simu¬ 
lation.  An  equivalent  non-rolling  linear  planar  simulation  was  run  and  the  nonlinear  simulations  were  run 
with  the  same  aerodynamics  and  autopilot  with  the  roll  rate  set  to  zero.  The  results  were  as  expected.  The 
nonlinear  model  reduces  to  the  linear  model  when  the  roll  rate  is  zero.  The  acceleration  and  fin  deflections 
for  the  three  loop  system  are  presented  in  Figures  3  and  4.  the  two  loop  responses  are  presented  in  Figures 
5  and  6 

A.  Roll  Rate  0.2  Hz  Commanded  Acceleration  0.5  G’s 

The  roll  rate  was  set  to  0.2  Hz  and  the  acceleration  command  was  set  to  0.5  G’s  in  the  inertial  z  direction. 
The  acceleration  time  histories  are  shown  in  Figure  7.  This  figure  shows  a  slightly  different  view  of  the  step 
response.  The  view  is  in  an  inertially  fixed  frame  that  is  initially  aligned  with  the  body  frame.  The  command 
is  in  the  inertial  z  direction  but  the  response  can  be  out  of  plane.  The  desired  response  is  to  remain  in  the 
inertial  x-z  plane.  As  can  be  seen  the  three  loop  autopilot  that  is  implemented  in  the  body  axis  system  is 
the  only  system  that  produces  significant  deviation  from  the  desired  response.  At  this  low  frequency  the 
deviations  might  be  deemed  acceptable  but  both  two  loop  and  the  roll  stabilized  three  loop  provides  a  better 
response. 

The  area  at  the  end  of  the  maneuver  is  expanded  and  is  shown  in  Figure  8  and  expanded  even  further 
in  Figure  9.  The  three  loop  implemented  in  the  body  axis  response  also  results  in  a  steady  state  limit  cycle 
of  approximately  the  same  magnitude  as  the  three  loop  roll  stabilized  response.  Notice  that  there  is  a  small 
error  in  the  two  loop  responses. 
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Figure  4.  Roll  Stabilized  Pitch  Fin  Deflection 


1.4 


1.2 


1 

g 

c  0.8 
o 

ra 

S 

8  0.6 
(J 
< 

0.4 


0.2 


0 

0  1  2  3  4  5 

Time(sec) 


Roll  Stabilized  Acceleration 


- 

/ 

- Body  Two  Loop’Roll  Stabilized  Two  Loop 

-  -  Non-rolling  Two  Loop 

Figure  5.  Roll  Stabilized  Acceleration 
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Roll  Stabilized  Pitch  Fin  Deflection 


Figure  6.  Roll  Stabilized  Pitch  Fin  Deflection 
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Figure  7.  Achieved  Acceleration 
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Figure  8.  Expanded  Achieved  Acceleration 
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Figure  9.  Further  Expanded  Acceleration 
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Roll  Stabilized  Pitch  Fin  Deflection 


Figure  10.  Roll  Stabilized  Pitch  Fin  Deflection 


The  roll  stabilized  fin  deflections  in  the  primary  channel  for  both  three  loop  autopilots  are  compared 
with  the  non  rolling  deflection  that  would  produce  the  planar  acceleration  are  shown  in  Figure  10. 

The  acceleration  responses  in  the  primary  direction  for  the  three  loop  autopilots  is  shown  in  Figure  11. 
The  two  loop  designs  are  shown  in  Figures  12  and  13.  The  off  axis  fin  deflections  are  presented  in  Figure  14. 


B.  Equivalence  of  Two  Loop  Closed  Loop  System 


As  can  be  seen  by  the  results  the  two  loop  autopilot  topologies  are  exactly  the  same  because  of  the  symmetry 
of  this  example.  For  the  two  loop  autopilot  based  on  the  body  axis,  the  fin  deflections  (removing  the 
measurement  notation)  are 

SY  =  KAy  (Ay  -  KssAVc)  +  KrR 
8p  =  KAz  ( KssAZc  —  Az)  +  KqQ 

and  for  the  roll  stabilized  system 


with 


and 


8y  —  A Ay  ^ Ay  KssAy\  R-  I\rr 

8p  =  Kaz  ( KssAZa  -  AAj  +  Kqq 

5y  =  Sy  cos(</>)  +  5p  sin(</>) 

Sp  =  —5y  sin(</>)  +  Sp  cos(</>) 

Ay  =  Ay  cos (cf))  +  Az  sin((/)) 

A,  =  —Ay  sin(</>)  +  Az  cos(^) 

q  =  Qcos((/))  +  Rsin((f>) 
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Figure  11.  Roll  Stabilized  Acceleration 
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Figure  12.  Roll  Stabilized  Fill  Deflection 
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Roll  Stabilized  Acceleration 


Figure  13.  Roll  Stabilized  Acceleration 
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Figure  14.  Roll  Stabilized  Yaw  Fin  Delflection 
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Achieved  Acceleration 


Figure  15.  Achieved  Aceleration 


r  =  — Qsin(^>)  +  R  cos(</>) 

Going  through  the  algebra  it  can  be  shown  the  two  loop  designs  are  exactly  the  same.  The  three  loop  designs 
are  not  because  of  the  integrator  in  the  loop. 

C.  Roll  Rate  2.0  Hz  Commanded  Acceleration  0.5  G’s 

The  roll  rate  was  set  to  2.0  Hz  and  the  acceleration  command  was  set  to  0.5  G’s  in  the  inertial  z  direction. 
The  acceleration  time  histories  are  shown  in  Figure  15.  As  can  be  seen  the  three  loop  autopilot  implemented 
in  the  body  axis  system  is  significantly  worse  than  the  other  three  autopilots. 

The  area  at  the  end  of  the  maneuver  is  expanded  and  is  shown  in  Figure  16.  Notice  how  the  axis  are 
of  a  different  scale  then  the  previous  plot.  The  limit  cycle  magnitude  has  increased  for  all  the  autopilots, 
including  the  body  based  three  loop  autopilot  which  is  not  shown. The  same  set  of  results  as  previously 
presented  are  shown  in  Figures  17  to  21. 

D.  Roll  Rate  2.0  Hz  Commanded  Acceleration  5.0  G’s 

It  was  desired  to  see  the  effect  of  increasing  the  commanded  G’s.  The  roll  rate  was  set  to  2.0  Hz  and  the 
acceleration  command  was  set  to  5.0  G’s  in  the  inertial  z  direction.  The  acceleration  time  histories  are  shown 
in  Figure  22.  As  can  be  seen  the  overall  response  is  slightly  different. 

The  area  at  the  end  of  the  maneuver  is  expanded  and  is  shown  in  Figure  23 

E.  Roll  Rate  2.0  Hz  Non-Zero  X  Moment  of  Inertia 

It  was  desired  to  investigate  the  effects  of  having  non  zero  x  moment  of  inertia.  The  nonlinear  equations 
were  used  instead  of  the  simplified  equations.  The  results  are  presented  in  Figures  24  through  26. 
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Figure  16.  Acieved  Acceleration 
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Figure  17.  Roll  Stabilized  Pitch  Fin  Deflection 
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Figure  18.  Roll  Stabilized  Acceleration 
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Figure  19.  Roll  Stabilized  Pitch  Fin  Deflection 
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Figure  20.  Roll  Stabilized  Acceleration 
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Figure  21.  Roll  Stabilized  Yaw  Fin  Deflection 
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Achieved  Acceleration 
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Figure  22.  Achieved  Acceleration 
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Figure  23.  Expanded  Achieved  Acceleration 
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Figure  24.  Achieved  Acceleration 
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Figure  25.  Expanded  Achieved  Acceleration 
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Figure  26.  Further  Expanded  Achieved  Acceleration 


F.  Very  High  Roll  Rates 

The  effect  of  increasing  the  roll  rate  significantly  beyond  the  bandwidth  of  the  closed  loop  system  will  now 
be  investigated.  Fist  we  must  determine  the  closed  loop  transfer  functions  of  the  non-rolling  airframe  and 
controller.  The  closed  loop  transfer  functions  can  be  determined  from  the  open  loop  system  and  the  feedback 
gains.  The  two  loop  acceleration  achieved  to  acceleration  command  is  given  by 

A  _  151.2013 

~AC  ~  s2  +  15.8891s +  151.2013 

151.2013 

(s  +  7.9446  +  9.3854*)  (s  +  7.9446  -  9.3854*) 

Since  this  is  a  second  order  speed  of  response  is  determined  by  the  distance  of  the  complex  poles  from  the 
origin.  The  closed  loop  system  bandwidth  is  approximately  2.0  Hz.  The  three  loop  systems  closed  loop 
transfer  function 

A  _  1054.5 

~AC  ~  s3  + 18.7s2  +  200.3s +  1054.5 

1054.5 

(s  +  4.6814  +  9.5242i)  (s  +  4.6814  +  9.5242*)  (s  +  9.3623) 

The  bandwidth  of  this  system  is  approximately  1.8  Hz.  The  roll  frequency  is  increased  to  10  Hz.  The  results 
are  somewhat  surprising  in  that  both  three  loop  topologies  degrade  significantly.  The  two  loop  has  degraded 
but  not  nearly  as  much  as  the  three  loop.  The  x  moment  of  inertia  is  again  neglected  and  the  commanded 
accelerations  returned  to  0.5  G’s.  The  acceleration  results  are  presented  in  Figure  27.  The  x  moment  of 
inertia  is  again  set  to  5  and  the  results  are  show  in  Figures  28.  These  show  the  three  loop  body  based 
autopilot  is  not  performing,  the  higher  roll  rate  just  makes  the  response  further  from  the  desired.  The  roll 
stabilized  three  loop  autopilot  is  degrading  to  an  almost  neutrally  stable  solution  by  10  Hz  and  has  crossed 
over  to  unstable  by  20Hz.  The  two  loop  autopilot  has  significantly  better  performance  at  these  higher  roll 
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Figure  27.  Achieved  Acceleration 


rates  with  the  x  moment  of  inertia  present,  even  though  it  does  not  produce  a  limit  cycle  around  the  desired 
set  point. 


X.  Conclusion 

Four  autopilot  topologies  were  examined  for  control  of  a  rolling  airframe.  These  were  the  classical  three 
loop  implemented  in  the  body  axis  system,  the  classical  three  loop  implemented  in  a  roll  stabilized  system 
and  a  two  loop  autopilot  implemented  in  both  the  body  and  roll  stabilized  frames.  It  is  shown  the  two 
loop  implementations  are  identical.  There  is  no  advantage  of  using  the  roll  stabilized  coordinate  frame  for 
the  two  loop  design.  The  three  loop  autopilots  produce  significantly  different  results.  For  this  example, 
all  autopilot  formulations  provided  adequate  responses  for  a  slowly  rolling  airframe.  The  example  had  the 
body  roll  rate  at  one  tenth  the  bandwidth  of  the  closed  loop  system.  The  body  based  autopilot,  even  at  this 
low  roll  rate,  showed  deviation  from  the  desired  response.  The  deviation  of  this  implementation  continued 
as  the  roll  frequency  increased  .  Under  no  conditions  did  the  system  become  unstable.  The  roll  stabilized 
three  loop  system  performs  well  while  the  roll  frequency  remains  below  the  closed  loop  bandwidth  of  the 
airframe.  At  ultra-high  (ten  times  the  closed  loop  bandwidth)  roll  rates  the  system  went  unstable.  The  two 
loop  autopilot  performed  better  at  the  ultra  high  frequencies.  This  system  did  not  show  signs  of  instability 
until  the  airframe  roll  rate  was  greater  than  50  Hz.  Zero  steady  state  error  was  not  achieved  for  the  two  loop 
autopilot  but  it  remained  small  under  all  conditions  except  when  the  airframes  roll  rate  was  larger  than  50 
Hz.  It  should  be  noted  a  first  order  160  Hz  control  actuation  system  was  included  in  all  models  in  order  to 
remove  the  algebraic  loop  in  the  simulation. 
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Figure  28.  Achieved  Acceleration 
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